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Abstract 

We propose that relative variability on short time-scales of the multiple images 
of a lensed quasar, after removal of the time delay, may be caused by hot spots 
or other moving structures in the accretion disk crossing microlens caustics caused 



'. by stellar mass objects in the lensing galaxy. Such variability has been reported 



in the two images of 0957+561. The short durations would be due to the high 



' rotation speed of the disk (v/c ~ 0.1), rather than planetary mass objects in 

I ' the slowly moving {v/c ~ 10" ) lens. This interpretation could be confirmed 

■ by finding periodicity, or correlations of the spectral and flux variations due to 

■ the Doppler effect in the disk. We also propose another signature of stationary 
accretion disks (with no intrinsic variability): the gradient of the magnification over 
the accretion disk should cause a relative color change between the images whose 

^ i sign and amplitude are correlated with the time derivative of the flux difference 

I between the images. Other color terms induced by the radial variation of disk 

colors are of second order in the magnification gradient. The methods proposed 
here can be used first to verify that accretion disks near supermassive black holes 
are the source of the continuum radiation from quasars, and then to study them. 
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1. Introduction 



For the past 30 years, the dominant view has been that the bulk of the lumi- 
nosity of quasars originates in accretion disks spiraling in towards super-massive 
black holes (Lynden-Bell 1969; see Rees 1984, Blandford et al. 1990, Lin & Pa- 
paloizou 1996 for reviews). It is therefore striking that to date, there has been no 
clear observational proof of this hypothesis. The main features of the theory are: 
1) The mass of the black hole is roughly estimated from its Eddington luminosity: 
M ^ 1O^M0(L/1O'*^L0). 2) The temperature of the accretion disk is estimated 
from the "blue bump" in the continuum spectrum to be T ~ 5 x lO'^K. 3) The 
characteristic radius of the accretion disk is estimated assuming a roughly black 
body intensity and is ~ 10^ AU for a bright quasar. The main observational dif- 
ficulty is that the angular size of the accretion disk at a cosmological distance is 
then very small (~ 1 //as), and so has not been possible to resolve to date. 

Microlensing provides a potentially powerful probe of quasars since the Einstein 
ring of a stellar-mass lens is typically ~ 10"^ AU, i.e. of the same order as the 
expected size of the accretion disk. Microlensing should occur in most lensed 
quasars. When multiple images of a quasar are produced by a lensing galaxy (with 
typical separations of 1"), the total surface density of the lens must be near the 
critical surface density. If a large fraction of this surface density is contributed 
by stars, then these stars will no longer cause microlensing lightcurves as isolated 
point masses, but instead their caustic curves are linked in a complicated network. 
The lightcurve of a source moving behind a lens then results from the variation 
of the magnification of several "microimages" with separations of the order of the 
Einstein radii of the stars, ~ 1 //as (e.g., Wambsganss, Schneider, & Paczyhski 
1990), and is therefore very complex. Every time the source moves inside (outside) 
a caustic two new images appear (disappear), and the magnification reaches a 
maximum that depends on the size of the emitting region in the source. 

Microlensing was first detected by Irwin et al. (1989) in Huchra's Lens (QSO 
2237+0305; Huchra et al. 1985), and shortly thereafter by Schild & Smith (1990) 
in the original double quasar, 0957-1-561 (Walsh et al. 1979). To verify that photo- 
metric variation is due to microlensing, it must be distinguished from other sources 
of variation, intrinsic variation in particular. If one of the macroimages is found 
to vary relative to the other, then one may conclude that one or the other image 
(or both) are being microlensed. However, the two images arrive at the observer 
with a relative time delay. Thus, the flux ratio of the two images can be measured 
only by comparing observations that are separated by this delay and, of course, 
this can only be done after the time delay has been measured. 

If the timescale of a microlensing event is long compared to the time delay, the 
latter can be ignored to first approximation. This is the case for 2237-1-0305: the 
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time delays among the four images are known to be of order days from the model of 
the lens, even though they have not been measured. Since the microlensing events 
are of order months, the time delay is not important. Schild & Choflin (1986) first 
measured the time delay for 0957+561 at r = 1.1 yr. Schild & Smith (1990) used 
this delay when they showed that the B image had brightened relative to the A 
image by about 0.2 mag over 10 years. This rise has subsequently slowed to a stop 
(R. Schild 1996, private communication). 

Intensive monitoring programs were then carried out to measure the time-delay 
more precisely. Since the quasar varies on scales of one or two days at the few 
percent level, a precision measurement did not at first appear difficult. However, 
it was found by Schild (1990) that there does not exist any value for r that would 
cause the 1-day to 3-month variations in the A and B images to coincide. He 
and his collaborators concluded that the underlying cause of most of these short 
timescale, low amplitude (Am ~ 0.02) variations is microlensing, not intrinsic 
variability (Schild 1990; Schild & Smith 1991; Schild & Thomson 1995; Schild 
1996). Kundic et al. (1996) also reported finding such variations in the course of 
measuring a more accurate time delay, r = 1.14 ± 0.01 yr. 

The microlensing interpretation of the short timescale events proposed by 
Schild has not received wide acceptance principally because the short timescales 
would seem to imply that the dark matter in the lensing galaxy is almost entirely 
composed of "rogue planets". The argument for the mass scale is simple. The 
characteristic physical Einstein radius is 



where i^oL; Dqq, and Dls ^'I'^ the angular diameter distances between the observer, 
lens, and source. For M — IMq, and adopting Hq — 55kms~^Mpc~^ and — 1, 
rE/DoL = (6.2, 2.1) //as for the lenses 2237+0305 and 0957+561, respectively. The 
crossing time is ts — te/vii where vt is the transverse velocity: 



Here, Vq, v^, and v^. arc the transverse velocities of the observer, lens, and source, 
and zi and Zg are the rcdshifts of the lens and source (Kayser, Rcfsdal, & StabcU 
1986; this formula is valid for flat cosmological models, in general the observer's 
velocity relative to the CMB frame must be multiphed by a factor that depends on 
the space curvature, see Miralda-Escude 1991). For 2237+0305 these parameters 



are (z/, Zg) = (0.039, 1.69) and (Dql, ^OS, ^ls) = (0.037, 0.29, 0.28)c//o"\ and for 




(1.1) 
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0957+561 they are {zi, Zg) = (0.36, 1.41) and (L>oL, ^os, ^Ls) = (0.21, 0.30, 0.18)cHq^. 
The Sun's motion relative to the cosmic microwave background is Vq — 300kms~"^, 
and if the transverse velocities of the lens and the source are similar then the 
timcscale to cross the Einstein radius for M = IMq is tE — (10, 30) yr for 
2237+0305 and 0957+561, respectively. Thus, when Irwin et al. (1989) reported a 
microlensing event with a duration of 3 months in 2237+0305, the naive interpre- 
tation was that this had been caused by a planet with a Jupiter mass. For events 
with durations of a few days in 0957+561, Earth masses would be inferred for the 
lensing objects. Wambsganss et al. (1990) pointed out that when the microlensing 
optical depth is large, events with timescales much smaller than tE are possible 
(caused, for example, by the passage of the source near a cusp catastrophe), so low- 
mass lenses are not automatically implied from the observation of just one event. 
However, the average density of caustic crossings should still not be larger than 
~ 1 per Einstein radius, so events on short timescales should be very rare unless 
planetary mass objects account for a large fraction of the lens surface density. 

However, an additional difficulty with the planetary-mass interpretation of the 
short events is the small size implied for the source. With the expected transverse 
speeds given above, the quasar in 2237+0305 moves by only 0.2 /las relative to the 
lensing galaxy over three months, corresponding to ~ 300 AU at the redshift of 
the quasar. This small size already implies a rather high temperature for the disk 
emitting the continuum {T ~ 10^ K; Ranch & Blandford 1991, see also Jaroszynski, 
Wambsganss, & Paczyhski 1992). For the shorter events in 0957+561, and with an 
expected angular velocity ~ 10 times smaller than in 2237+0305 (due to the larger 
distance to the lens), the inferred size of the source would be highly implausible 
within the framework of current models. 

In this paper, we propose that microlensing can be used to resolve the structure 
of accretion disks. In § 2 we examine the possibility that short timescalc events, 
of the type reported by Schild (1996), could be produced by intrinsic variability 
in the accretion disk combined with the effects of microlensing. In § 3 we turn to 
the case of a disk with no intrinsic variability, and propose a method to demon- 
strate the existence of an accretion disk from the spectral changes in the different 
macroimages as microlensing events take place. 
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2. Short Timescale Variability 



The main difficulty in explaining the short timescale events is the small proper 
motion expected, corresponding to the typical galaxy velocities generated by large- 
scale structure of ~ 300 kms~^. However, the orbital velocity of the accretion disk 
in the continuum emitting region should be much larger. For example, if the black 
hole mass is M = 10® M©, the circular velocity is 3 x lO^kms^^ at a radius of 
100 AU, implying that the proper motion of any orbiting blob is (5,40) times 
larger than the proper motion of the lensing galaxy in 2237+0305 and 0957+561 
[from cq. (1.2)]. We therefore propose that the short events observed in 0957+561 
(and possibly some of the variations in 2237+0305 as well) are caused by "spots" 
moving at the typical orbital velocities. Such spots might be caused by a variety 
of phenomena, including instabilities in the accretion disk that create "hot spots" 
which are in roughly circular orbits, gas clouds or stars that fall to and crash 
against the accretion disk, or blobs ejected by a jet near the black hole. They 
could even be "cold spots", relatively confined subluminous regions like sunspots 
on the Sun. 

If the magnification of each image were uniform over the whole accretion disk, 
then the intrinsic variability of the quasar caused by such spots would be the same 
in all the images after correction for the time delay, and therefore they would not 
be assigned to microlensing. But if the magnification varies over the region of the 
motion of a hot spot by AA due to microlensing, and the hot spot contributes 
a fraction AF to the total quasar fiux F, then we expect a relative variation of 
(AA/A){AF/F) of the flux in different images. The timescale of these short events 
should be shorter than or of order the dynamical time in the accretion disk; for a 
radius r = 100 AU and velocity f /c ~ 0.1, this is about one month (which should 
be redshifted by a factor ~ 2.5 for the two quasars we have mentioned). If the 
spot is a transient phenomenon (changing its fiux AF over a dynamical time) , we 
should expect an intrinsic variation of AF/F, so the intrinsic and relative variation 
should be correlated. 

The largest variations in magnification will occur when a caustic transits the 
accretion disk. Two images of the accretion disk will then be merging on a criti- 
cal line and will be highly magnified. Let us assume that the spot has a surface 
brightness higher by a factor (1 + fg) compared to the average of the accretion 
disk {fs > for hot spots and fs<0 for cold spots). Then the radius of the 
spot is Rg = (AF/(/sF))^^^ r, where r is the radius of the disk. Since the maxi- 
mum magnification on a fold catastrophe is proportional to Rg (e.g., Schneider, 
Ehlers, Sz Falco 1992), the fractional contribution of a spot to the total fiux can 
rise to {AF/ F){r/ Rg)^/^ = (| AF|/F)3/4 V4. por example, a spot from a region 
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contributing only 0.2% of the total light and with ~ 1 could cause fluctuations 
~ 1%. 

For a spot that is either falling to or being ejected from the black hole, the 
caustic will be crossed once. But if a spot is in circular orbit and can survive 
for several orbits, then the caustic will be crossed repeatedly. An example of the 
lightcurve that may result from this is illustrated in Figure 1. The middle panel 
shows the lightcurve for a single hot spot in circular orbit. We assume that the 
caustic crossing the disk is on a sufficiently large scale, so that the magnification 
over the disk can be approximated as ^4 = + K/x^/'^. There is a characteristic 
pattern of variations as the spot crosses the caustic: an "M-shaped" event is pro- 
duced at every orbit of the spot. These events might therefore be confused with 
similar events when the motion is linear and two folds of the same closed caustic 
are crossed. The periodic events become shorter until they disappear (the time- 
reversed lightcurve is of course equally likely if the motion of the disk is reversed). 
Such a spot would be likely to have a different spectrum than the average for the 
quasar (it might, for example, be much closer to a blackbody). The spectrum is 
redshifted and blueshifted as the spot moves along the orbit (the radial velocity 
is shown in the lower panel for our example). There should then be a correla- 
tion between the spectrum and the magnification. In our example, the spot is 
always more blueshifted at the beginning of each event. This blueshift effect may 
be detectable even in the continuum spectrum of a spot, given the large velocities 
involved. Notice that such changes of the spectrum due to microlensing of spots 
can be separated from any other intrinsic changes of the spectrum of the entire 
disk, if one focuses on relative variations of the spectrum between the multiple 
images of the quasar after correcting for the time-delay. 

Whether an orbiting spot can survive for several orbits depends on the process 
that generates it. A topologically confined spot, such as a cold spot due to a 
magnetic vortex, may persist for many orbital periods. Long lived spots may also 
be generated by gravitationally confined disturbances in the gaseous disk, and by 
the effects of compact objects orbiting inside the disk (such as black holes or other 
massive stars that may have fallen to the disk). On the other hand, hydrodynamical 
disturbances would expand in their sound crossing time and would then be sheared 
out into arcs by differential rotation. For example, if the ratio of the sound speed 
to the circular velocity is 0.002 (e.g. T = 5 x 10^ K and v/c = 0.1), a spot with 
orbital radius r should, after periods, have expanded radially to a size 0.01A"r 
due to expansion, and tangentially to a size O.lA^^r due to differential rotation. 
Therefore, such a spot cannot survive for more than a few periods. 

In brief, spots can give rise to a variety of short timescale signals in the differ- 
ence between the flux from two quasar images. The deviations in the microlensed 
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image minus the non-microlensed images can be both positive and negative and 
can contain repeating, non-repeating, and quasi-repeating signals. 

3. Long Timescale Variability from a Steady Disk 

The color of a microlensed accretion disk will in general be different from that 
of the disk in the absence of microlensing, due to differential magnification over 
the disk. If there are two or more macroimages of the disk, then this color shift 
can be determined by comparing the colors of the two images observed at two 
times separated by the measured time delay. The effect of microlensing can then 
be unambiguously distinguished from intrinsic color changes. Color differences 
between the images will arise if either or both are microlensed but for simplicity 
of discussion we will assume that one image is microlensed and the other is not. 

Color changes are caused by two distinct effects of different orders. Both 
require that the magnification vary over the source. The variation can be either 
smooth or due to a caustic. We focus here on the case of smooth variations in the 
magnification because these are present generically. 

The first order effect arises from a couphng of the hnear color gradient across 
the accretion disk (due to the disk's rotation) with the gradient of the magnifi- 
cation. This effect is analogous to the line shift induced by microlensing of stars 
which has already been studied in some detail (Maoz & Gould 1994). In the stellar 
case, the stellar lines are broadened by rotation because one side of the star is 
rotating toward the observer and is blueshifted while the other side is redshifted. 
If, for example, the blueshifted side of the star is more highly magnified, the blue 
wing of the line will also be more highly magnified and the line centroid will shift 
toward the blue. For a star with projected rotation speed vsini and radius r, the 
shift is 

( 

Ai> = -r|Vln A|vsinisin7, (3.1) 

where V In A is the logarithmic gradient of the magnification and 7 is the angle 
between the magnification gradient and the projected axis of rotation. The factor 
C depends on the details of the geometry. For a line emitted uniformly over the 
surface of a star, C = 1, while for a ring of emission such as the Ca II line at 393 
nm in giant stars (Loeb & Sasselov 1995), C = 2. 

The origin of the broad lines seen in virtually all quasars is a subject of consid- 
erable debate. From their width (v ~ 10^ kms^^), they must come from regions 
that are at least ~ GM/v^ ~ 10^ AU from the center. If they come from the ac- 
cretion disk (or from any other structure with organized motion on similar scales) 
then they will be subject to a line shift given by equation (3.1). While the exact 
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value of the shift will depend on the details of the geometry and will in addition 
be time dependent, the general order of the effect should be ~ i'/4 ~ 2000 kms~^, 
assuming r|Vln74| ~ 1. Notice that for r = 10^ AU, the angular size is ~ 1 //as; 
since the typical Einstein radius is only 1 to 5 times larger, a magnification gradi- 
ent |VlnA| ~ 1/r should be common. Hence, there should be an observable shift 
at some times if the broad line region is indeed associated with the accretion disk. 
Schild & Smith (1991) measured the Mgll line of the two images 0957+561 in two 
observations separated approximately by the time delay. No difference in the lines 
was detected although it is not clear that the data were of sufficient quahty to see 
the predicted line shift. In any event, a definitive test of the accretion-disk origin 
of the broad lines would require measurements at multiple epochs. 

Even if the accretion disk does not give rise to broad lines, a magnification 
gradient can still generate a color shift. In this case, the entire quasi-thermal 
spectrum ( "blue bump" ) may be loosely considered as a giant "emission line" with 
a fractional width of order unity. More concretely, consider a color formed from 
the ratio of fluxes at two broad-band wavelengths and Xr with spectral slopes 
as and a^. By a calculation analogous to the one of Maoz & Gould (1994), the 
first order color shift A{B — R)i is given by 

2.5 ,, sinisin7 (rufoiR — a??)) ,„ ^, 

A{B - i?)i = ^ iVln^l ^ (3-2) 

where {rv{aB — o:r)) is the intensity weighted mean over the disk profile, and 7 
is the angle between the minor axis of the projected disk and the gradient of the 
magnification. Thus, for \aB — oir] ~ 1, VlnA ~ 10~^AU~-^, and characteristic 
disk radius ~ 300 AU, one might typically expect color shifts of order ~ 1%. 

The second effect that causes color variations in microlensing on a moving 
stationary disk, which is of second order, is due to the radial color variation. 
This was investigated numerically by Wambsganss & Paczyhski (1991). The color 
change can be written, to second order, as (see, e.g., Gould & Welch 1996) 

A{B - R)2 ~ ^^^-^{rVlnAf, (3.3) 
o 

where r is the characteristic radius of the disk, and 



2 /Q°°rfr'r'3^B(rO 
Jg°° dr' r' Sb{t') 



As = 0^° 00 J , (3.4) 



is the second radial moment of the disk in the B band (with the same definition 
for the it! band). If, for example, A^ — ~ 0.1, and (rVln^l) ~ 0.3 (which we 
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would expect typically), this color term is still substantially smaller than the first 
order term given by equation (3.2). The first order term, which is proportional 
to V In A, then dominates the relative color variations of the quasar images. The 
time derivative of the magnification of the microlensed image is also proportional 
to V In A: 

dlnA Dos COS0 

—r—^vt\V\nA— — — — , (3.5 

at Dql ^ + zi 

where is the angle between the magnification gradient and the transverse veloc- 
ity. Hence the first order color term should be directly proportional to the time 
derivative of the magnification ratio (after taking out the time delay). This allows 
it to be distinguished from the second order term. In practice, the proportionality 
between the color and the magnification variation should not be perfect because, 
as the magnification varies, the angle should also vary. Moreover, some of the 
color variations may come from a sufficiently large region in the disk so that the 
first-order term does not dominate. The maximum scale at which a first-order 
expansion of the magnification can be applied is reduced, of course, if low-mass 
stars or brown dwarfs are highly abundant, and is also highly variable in the 
magnification patterns generated by the random superposition of stars. Neverthe- 
less, a correlation between the color and magnification variation should clearly be 
present. Detailed numerical simulations of lightcurvcs, with realistic disk models 
around Kerr black holes (e.g., Jaroszynski et al. 1992) should help determine the 
signatures of accretion disks that can be found. 
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FIGURE CAPTIONS 



1) A representative example of the lightcurve that may result from an orbiting 
hot spot. The upper panel shows the trajectory of the spot over one orbit. 
The thick line represents the caustic (which is moving upwards relative to the 
disk) and the dashed line is the line of nodes of the orbit. The magnification 
near the caustic is assumed to be ^4 = 1 + (&/?/) where y is the distance 
to the caustic, and we choose the constant b — 2 //as, equal to the Einstein 
radius of a 1 Mq star. The radius of the orbit is r = 100 AU, inclined at 
an angle 50° to the line-of-sight, and the line of nodes has an inclination of 
40° relative to the caustic. The orbital velocity of the spot is assumed to 
be v/c = 0.1, and the relative proper motion of the caustic and the disk 
is 40 times smaller than that of the hot spot. The solid line in the middle 
panel shows the magnification of the spot as a function of time, for a size 
of the spot equal to 0.1 times the radius of the orbit (this size determines 
the maximum magnification at each caustic crossing). The magnification 
can be measured from the observed flux variation relative to another, non- 
microlensed image after correction for the time-delay. The dotted line is the 
actual flux variation of the spot, which is affected by intrinsic variation due 
to Doppler effects when the spot is assumed to emit a constant flux in its 
rest-frame. The lower panel shows the radial velocity of the spot. 
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